import sparse
from sparse._compressed import CSR, CSC

gcxs_a = sparse.random((10000, 20000), .01, format="gcxs")
gcxs_b = sparse.random((10000, 20000), .01, format="gcxs")

csr_a = CSR(gcxs_a)
csr_b = CSR(gcxs_b)

scipy_a = csr_a.to_scipy_sparse()
scipy_b = csr_b.to_scipy_sparse()

%timeit gcxs_a + gcxs_b
# 999 ms ± 10.2 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)

%timeit csr_a + csr_b
# 32.2 ms ± 1.9 ms per loop (mean ± std. dev. of 7 runs, 10 loops each)

%timeit scipy_a + scipy_b
# 27.5 ms ± 894 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)